Response to perturbations for granular flow in a hopper 
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We experimentally investigate the response to perturbations of circular symmetry for dense gran- 
ular flow inside a three-dimensional right-conical hopper. These experiments consist of particle 
tracking velocimetry for the flow at the outer boundary of the hopper. We are able to test com- 
monly used constitutive relations and observe new granular flow phenomena that we can model 
numerically. Unperturbed conical hopper flow has been described as a radial velocity field with 
no azimuthal component. Guided by numerical models based upon continuum descriptions, we 
find experimental evidence for secondary, azimuthal circulation in response to perturbation of the 
symmetry with respect to gravity by tilting. For small perturbations we can discriminate between 
constitutive relations, based upon the agreement between the numerical predictions they produce 
and our experimental results. We find that the secondary circulation can be suppressed as wall 
friction is varied, also in agreement with numerical predictions. For large tilt angles we observe the 
abrupt onset of circulation for parameters where circulation was previously suppressed. Finally, we 
observe that for large tilt angles the fluctuations in velocity grow, independently of the onset of 
circulation. 

PACS numbers: 47.57.Gc,45.70.-n,83.80.Fg,47.85.M- 



I. INTRODUCTION 

Many practical processes involve the flow of dense 
granular materials. Applications range from food grains, 
ores, and coal to pharmaceutical powders. To predict 
flows of dense granular states, one would like a homog- 
enized/continuum model that reflects the complexity of 
grain-scale interactions but can be applied at much larger 
scales. Continuum approaches frequently describe gran- 
ular materials using elasto-plastic models in which the 
material can support shear stresses up to a certain value 
before yielding and deforming irreversibly. In the sim- 
plest case, this approach amounts to Coulomb's law of 
friction. 

The most common yield criterion is the Mohr-Coulomb 
criterion, which proposes a linear relationship between 
the magnitudes of normal stress a and shear stress r 
in two-dimensions. Because the choice of axes dictates 
the relative magnitude of the normal and shear stresses, 
the range of possible axes maps out a circle in the space 
of normal and shear stress — the "Mohr circle". The 
diameter of this circle is proportional to the magnitude 
of the stresses [![. 

In the normal-shear stress space, the Mohr-Coulomb 
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yield criterion is the straight line r = \io + c where \x 
is the coefficient of friction and c is the cohesion. It 
follows for cohesionless granular materials, the case we 
consider here, that the yield criterion must pass through 
the origin. The stress within the material, and hence 
the diameter of the Mohr circle, can increase only until 
the yield criterion is tangent to the circle. The point of 
contact then indicates an axis of a coordinate frame in 
which the ratio of shear to normal stress (equivalent to 
the mobilization of friction) is maximized. In this two- 
dimensional analysis of a three-dimensional material, the 
direction of the maximally mobilized shear stress and the 
third, neglected, axis span the "mobilized plane" along 
which motion of the granular material is most likely in 
response to stress [l|. 

There are several approaches that have taken for gen- 
eralizing the Mohr circle picture to three dimensions. An 
example is the Matsuoka-Nakai criterion. By geometri- 
cally averaging the three planes obtained by considering 
pairs of dimensions a "spatially mobilized plane" is ob- 
tained. This is similar to the generalization of the Tresca 
failure criterion for metals to three dimensions by von 
Mises [2j. Tri-axial experiments support the validity of 
the Matsuoka-Nakai criterion, indicating that it may cap- 
ture the average behavior of a granular material [3j. 

Unfortunately, the mathematics of continuum descrip- 
tions of granular materials are rife with instabilities, ill— 
posedness, and complex non-linearities [J]. This means 
that evaluating the response of granular materials to 
perturbations that cause unusual behavior in continuum 
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models is necessary to determine the physical relevance 
of such approaches. 

The right-conical hopper provides an ideal test system 
for granular flow in that it is not only an experimentally- 
accessible system, but also because the behavior is com- 
monly assumed to tend to a steady-state with known 
velocity and stress field solutions derived from soil me- 
chanics. The so-called Jenike radial solutions describe 
the velocity and stress fields in an infinite hopper as self- 
similar functions of radial position alone, without non- 
radial velocity components Q. In two dimensions, such 
soil mechanics descriptions of wedge hoppers have been 
studied extensively both experimentally (| 0, H, [§] and 
theoretically [l(| El- Three-dimensional, right-conical 
hopper flow has been studied through exp eriments and 
discrete element simulations in the past |12l. Il3j| . but these 
studies did not make direct comparison to soil mechanics 
predictions. 

Recent numerical work has taken advantage of the self- 
similar nature of the solution even in relatively general 
three-dimensional geometries such as "pyramidal" hop- 
pers. By perturbing the geometry of the simulated hop- 
per, the response of the continuum description can be 
examined. In particular, numerical work has found that 
secondary, non-radial circulation currents arise. These 
currents depend sharply upon model parameters in non- 
trivial ways OH Bit. 

We make quantitative experimental observations of 
a right-conical granular hopper to determine if these 
numerically-predicted behaviors, the result of non- 
linearities in the continuum description of granular ma- 
terials, have physical significance. We investigate both 
the influence of tilt-induced perturbations and changes in 
boundary roughness. By comparing with numerical pre- 
dictions we investigate the ability of different constitutive 
relations to correctly model experimentally observed phe- 
nomena. Since we obtain localized time-resolved veloc- 
ities, we also identify useful statistical characterizations 
of flowing dense granular matter. 



II. SOIL MECHANICS APPROACHES TO 
GRANULAR MATTER 

Before turning to the experiments, we give a brief sum- 
mary of the relevant soil mechanics constitutive relations. 
At issue are the determining equations for the stress ten- 
sor. Force balance gives six equations for the nine compo- 
nents of the three dimensional stress tensor. To close this 
system of equations we need additional constitutive re- 
lations describing the yielding, plasticity and flow of the 
material [Hj]. Consideration of the velocity field v de- 
scribing the flow introduces three more unknowns bring- 
ing the necessary number of additional constraint equa- 
tions to six. 

We use Levy's flow rule to introduce five constraints 
[l], [l3| ■ Missing is a final constitutive relation describing 
the stress sufficient to cause plastic rearrangement, or 



yield, within the material [la ]. 

In metal plasticity, the von Mises failure criterion is 
commonly used as a constitutive relation. The von Mises 
yield surface in stress space, however, is independent of 
the total pressure: 



(<Tl - er 2 ) 2 + (er 2 - er 3 ) 2 + (<T 3 - (Ti) z 



(1) 



where <j\, a%, 03 are the principal stresses denoted in 
order of magnitude and c is a constant. This is necessar- 
ily at odds with the Coulomb nature of granular failure. 
Replacing the constant term on the right-hand side by a 
term involving the principal stresses generalizes the von 
Mises criterion to granular materials by the introduction 
of a pressure-dependent yield surface: 

(<Tl - °2? + {V2 - T 3 ) 2 + (<7 3 - CT l) 2 = 6shx9 s 2 I* (2) 



or equivalently, 
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where the three stress invariants are I\ = 
^ (a 1 + 1T2 + 03)) (equivalent to the isotropic pres- 
sure), I 2 = - (0-20-3 + <73<Ti + 0102)) and 7 3 = o'lO'ijCTp 
and 9 S is the internal angle of friction of the sand [l|, [2| • 
This "granular von Mises" condition is commonly used 
as the missing granular constitutive relation. 

As an alternative, Matsuoka and Nakai [1J apply the 
Mohr-Coulomb failure criterion \<Jjk\ < t&n-4>i (J kk for a 
cohesionless material to obtain, for each pair of axes j 
and k, a Mohr circle with radii Ri = (<jj — cj.) sin(9 s and 
a yield criterion that makes an angle <pi at the origin such 
that cos 4>i = \<jj — (Tk\/Ri- They argue for the following 
constitutive relation: 
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which they show to be equivalent to: 

h * h/h = c 



(4) 
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Gremaud et al. have computed flows in relatively gen- 
eral geometries for both of the above plasticity models 
[HI . Ila ] . The authors model granular flow in a vertically 
oriented infinite right-conical hopper described in spher- 
ical coordinates where r is the distance from the origin, 
which is at the tip of the cone, 9 is the angle measured 
away from the axis of symmetry and <f> is the azimuthal 
position about the symmetry axis. The authors vary the 
wall angle azimuthally about the constant angle 9 W as 
9 = 9 W + ecosm^. The case m — 1 corresponds to elon- 
gating the cross-section of the hopper in a manner that 
is roughly similar to the effect of tilting the hopper with 
respect to gravity. Assuming e is small (between —5° and 
5° for a tilted hopper), the velocity can be expanded as 
v — v + ev 1 + ... giving the radial and azimuthal com- 
ponents of the velocity as: 

v r w u° + evl 

« ev\ (6) 



3 




FIG. 1: (Color online) Numerical simulations of the velocity 
stream functions for tilted hoppers predict the rise of two 
secondary circulation cells, rotating in opposite directions. 
Darker tones represent greater velocity. In the top half of 
the figure the circulation is clockwise while in the lower half 
the circulation is counterclockwise. The strength and shape 
of the cells depend on the angle of tilt, angle of internal fric- 
tion, and coefficient of wall friction. The azimuthal angle <f> 
is taken to be zero on the line of symmetry between the two 
cells. 



where an unperturbed (e = 0) hopper has Jenike radial 
flow. 

Two main observations were drawn from the numer- 
ical experiments of Gremaud et al. First, when verti- 
cal axi-symmetry is broken, secondary circulation takes 
place, i.e., becomes non-zero and leads to two counter- 
rotating cells perpendicular to the tilt axis, as in Fig. [T] 
Second, the computed circulation currents were much 
larger for a Matsuoka-Nakai material than for a von Mises 
one. 

The numerically predicted components of velocity 
and depend upon the choice of model parameters and 
constitutive relation. When using the Matsuoka-Nakai 
relation, v 1 ^ varies strongly as a function of wall friction, 
fi w , ranging from the same order as the radial flow to 
three or more orders of magnitude smaller. Similar be- 
havior is seen when the von Mises granular criterion is 
used, but the effect is smaller by orders of magnitude 

The numerical model assumes mass flow, where all 
the grains are moving and there are no stagnant regions 
within the hopper. For this reason, simulations were not 
conducted above fi w — tan(0 s ) = 0.5, where the angle of 
friction of the simulated material was 9 S = 30° [l4[ ■ 
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FIG. 2: Schematic of hopper in which we observe flowing sand. 
To perturb the flow we tilt the hopper up to 23 degrees. 



III. METHODOLOGY 

Our experimental test hopper is an approximately- 
conical, regular polyhedron composed of fourteen brass 
wedges (fi w = 0.34 when bare) and two clear Plexiglas 
wedges (p. = 0.51). The hopper, depicted in Fig. is 
~ 45 cm wide at the top and the walls are angled at 
6 W = 24°. Flow is slowed by a nozzle with a 0.6 cm 
diameter opening to ~ 6 g/s and the nozzle is screened 
from the hopper by a 3 cm diameter tube that is suffi- 
ciently long (~ 20 cm) for the influence of the nozzle to 
be neglected. Sand runs from the opening at the bot- 
tom of the hopper into the discharge tube, with foam 
wrapped around the joint to prevent spillage. It takes 
~ 20 minutes for the hopper to drain if not refilled. A 
second, feeder hopper is placed on axis at the top to 
provide refilling from an approximate point source ~ 14 
cm within the hopper. The feeder hopper is refilled by 
hand as needed and keeps the main hopper nearly-fixed 
at roughly one third full. By feeding from the small hop- 
per we maintain a well-controlled preparation of the sand 
in the primary hopper. 

We note that the 16-sided cross section of the hop- 
per involves some perturbation from circular symme- 
try. However, we estimate that these are 15 to 20 times 
smaller than the typical perturbation that we achieve by 
tilting the hopper. In addition, we expect that any small 
flows induced by these perturbations would not have a 
substantial swirling component. 

In these experiments we use Ottawa sand (ASTM C- 
190). The sand (coefficient of friction fi S and — 0.42, den- 
sity 1.64 ± 0.01 g/mL, and diameter 420 < d < 595 fim) 
is mixed with a small quantity of tracer sand that has 
been dyed dark red with instrument ink and repeatedly 
baked dry. 

The entire hopper can be physically tilted, perturbing 
the geometry with respect to gravity. The feeder hopper 
tilts with the system, staying on the axis of the main 
hopper. A digital (CCD) camera, bolted to the side of the 
hopper, tilts with the hopper and can be aligned with the 
Plexiglas "windows" in one of two configurations so that 
images showing the passage of sand may be recorded with 
a computer. The hopper is tilted so that the windows 
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are on the axis of the tilt, where v<j, is predicted to be 
maximal. A plumb line is hung within the hopper and 
imaged to determine the direction of gravity in the frame 
of the camera for each tilt. 

Since the hopper is made of flat wedges shaped into a 
cone, each wedge has a different orientation relative to 
the camera. To account for the tilt of the camera rela- 
tive to the orientation of the window wedges we mount 
the camera two different ways. In the first, we align the 
camera so that the tilts of the two windows are equal but 
opposite in orientation relative to the camera. In this ar- 
rangement, any geometric effect caused by the tilt of one 
window should be reversed in the other window. 

We use the second camera arrangement when we line 
the walls of the hopper with materials with different coef- 
ficients of friction. To minimize the impact of the viewing 
window when lining the hopper we cut a hole no bigger 
than the size of the field of view of the camera into the 
liner. In this case we align the camera flush with just one 
window of the hopper. 

The coefficient of friction, fi w , of a particular liner ma- 
terial is determined by using a weight and pulley system 
to drag a tray covered with the material across a bed of 
Ottawa sand. We capture the motion of the tray using a 
high-speed digital camera and analyze the images to de- 
termine the acceleration of the tray. For each run there 
is a period where the forces on the sled are roughly con- 
stant and we perform a linear fit to the velocity in this 
region to determine acceleration. By varying the load 
applied to the tray as well as the weight acting through 
the pulley we can account for systematic frictional forces 
and deduce the sliding coefficient of friction for sand on 
the liner material. 

Data sets for determining the velocity field typically 
consist of 1000 frames taken once every two seconds — 
covering roughly thirty-three minutes of flow. For a given 
combination of tilt and wall material, from two to six 
data sets are used to determine each plotted point for v. 

As illustrated by Fig. [31 particle tracking is performed 
by thresholding each frame into a binary intensity image 
where 'ones' identify pixels with intensity that is a tun- 
able number of standard deviations away from the local 
mean intensity in the region of the pixel. The thresh- 
olding is done by region to account for any non-uniform 
lighting. Clusters of 'ones' larger than a tunable size are 
then identified as the dark, tracer particles amidst the 
lighter sand, although occasionally bright or reflective 
sand grains are also resolved and tracked. Typically a 
dozen or more grains are tracked in a frame. 

Since the grains are slightly irregular, the cluster iden- 
tifying a given tracer grain changes as the particle ro- 
tates along the window surface. Typical tracer particles 
are roughly 20 x 20 pixels in size, though smaller dark or 
bright spots are tracked and used if possible. Each clus- 
ter is assigned a "center of mass" determined from the 
spatial distribution of its pixels, and then is identified 
with the nearest center of mass in the previous frame 
to construct a track of particle position in consecutive 



frames. If the nearest center of mass is outside a tunable 
maximum radius, the identification is rejected and a new 
particle track begins. Once a track is twenty-one frames 
long, a velocity is calculated using a least-squares fit and 
the interpolated velocity is assigned to the location of the 
tracer in the eleventh frame, as in Fig. 01 where out of 
the dozens of possible tracers only three are being used. 
If the velocity is zero, indicating an edge of the window 
is being tracked, it is ignored. Particle tracks begin and 
end as tracer particles move into and out of the field of 
view of the camera due to the overall flow and motion 
towards and away from the window — the mean track 
length is roughly 26 frames. In this manner we generate 
a time-integrated velocity field by binning velocities into 
regions by tracer location over the length of the run. 

Once a velocity field is generated for a particular com- 
bination of tilt angle and wall friction, as in Fig.O we can 
calculate the ratio v$/v r . If we have recorded using two, 
off-setting windows, we calculate the ratio separately for 
the two windows. We first linearly fit both components 
of velocity horizontally (fixed r) across the velocity field 
to the distance from the middle of the image to com- 
pensate for the slight tilt of the windows away from the 
camera. We interpolate using this fit to find the velocity 
components nearest the middle of the field of view of the 
camera from both windows, corresponding to <j) = ir/2 
(<p = is in the plane of the tilt). If we are using a single 
window, we perform a quadratic fit across the velocity 
field to account for (f> variation (a higher order effect in 
the two- window case) . We again interpolate the velocity 
components at <fi = tt/2. 

We find both the azimuthal and radial components of 
the velocity to vary radially as 1/r 2 , in agreement with 
mass conservation and the Jenike solutions. Given that 
both components vary as 1/r 2 , the ratio of the two ve- 
locity components, v$jv r , should be independent of r. 
Linear fits to 1/r 2 of the two velocity components are 
made and for each radial bin of the vector field the ratio 
of the v^-fit to the tv-fit is calculated and the mean is 
recorded as the observed ratio. Error bars are established 
using the standard deviation of these ratios. 



IV. RESULTS 

We find evidence for a non-radial component of the 
velocity field as a function of hopper tilt and wall friction. 
We also find that small perturbations of geometry, such 
as slight misalignment of the wedges forming our test 
hopper, causes a small amount of apparent secondary- 
circulation. With careful alignment and balancing of our 
untiltcd experimental hopper we were able to reduce, but 
not completely eliminate, this apparent non-radial flow 
even for the untilted case. In order to compensate for 
the observed small non-radial component of the flow, we 
rotate the experimental images to a frame in which the 
untilted hopper flow is perfectly radial. 

As the measurements from our test hopper shown in 





FIG. 3: Velocities are determined from identifying tracer particles and regressing over twenty-one consecutive frames to assign 
a velocity to the eleventh frame. Shown here are the actual frames (left) and thresholded tracer bitmaps (right) for a given 
frame (top) and the twentieth frame later (bottom), corresponding to 40 s of experimental observation. The dark dots indicate 
centers assigned to large tracer particles. 



Fig. [6] indicate, in a bare brass hopper (/j, = 0.34) sec- 
ondary circulation does in fact follow linearly with tilt 
angle for small perturbations (e < 0.15). At larger an- 
gles, the dependence appears to depart from linear, which 
might be expected, since higher-order terms should even- 
tually become important. Note that the ratio of az- 
imuthal to radial velocity, v^/v^t, follows systematically 
for tilts in both directions, indicating that our observa- 
tions are not an artifact. 

Numerically it is possible to predict v^/v^ as a func- 
tion of n scm d, 6 W aii and fi wa ii, although there are many 
domains where the predicted magnitude of secondary cir- 
culation is too small to be observed. For our system, 
wall friction is the most readily tunable parameter and 
we observed the flow for three additional wall frictions, 
as indicated by Fig. [7l 

By choosing the appropriate material, we were able 
to examine circulation in regimes where azimuthal flow 
would be suppressed. One such domain is that of low wall 



friction values, (such as Teflon — fi = 0.05). In the linear 
region (e < 0.1) there is no detectable secondary circula- 
tion, within experimental resolution, when the hopper is 
lined with teflon. As the tilt increases however, a strong 
azimuthal component is observed in the opposite direc- 
tion of the flow in the bare hopper. The tilt angle for 
the onset of flow showed some variation, possibly indi- 
cating hysteretic effects or some other influence beyond 
the controls of the experiment. 

As wall- friction is increased beyond HcriticaU soil me- 
chanics models of granular flow predict a transition from 
mass-flow — where all the material in the hopper is mov- 
ing — to funnel flow where there are non-moving, stag- 
nant regions at the wall. However, when we line the 
hopper with Emery paper, which has a higher coeffi- 
cient of friction fi = 0.75 than the interparticle friction 
fJ-crit = tan S and we continue to observe some radial flow 
at the window, although there is no azimuthal component 
within measurement error. 
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FIG. 4: The velocities assigned to the eleventh frame in 
the sequence depicted in Fig. [3] indicate that in that frame 
only three tracer particles had sufficiently-long unambiguous 
tracks to be used. 
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FIG. 5: A time- averaged velocity field is generated by spa- 
tially binning the velocities for every frame. 



Although extreme values of wall friction induce large 
changes in flow, we find that the flow is relatively robust 
to small variations in wall friction. Changing the wall 
friction slightly from the bare hopper (from fi w = 0.34 
to 0.26), by lining with low-friction Ultra-High Molecular 
Weight Polyethylene (UHMWPE) , does not qualitatively 
change the secondary circulation. 



FIG. 6: The ratio of azimuthal to radial velocity for the right- 
(•) and left-hand (o) windows for a bare (ji w — 0.34) hopper. 
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FIG. 7: Comparison of the ratio of azimuthal to radial velocity 
for Teflon {fj, w = 0.05 — O), UHMWPE (/j, w = 0.26 — □) 
and Emery paper (fj, w = 0.75 — A) to a linear regression of 
the bare (fi w = 0.34) hopper observations (solid line). 



V. COMPARISON WITH SOIL MECHANICS 
PREDICTIONS 

The dashed and dotted lines in Fig. [8] show azimuthal 
to radial velocities ratios as calculated for varying tilts in 
a bare hopper using a spectral method [l6| . We examine 
the ratio at the point in the hopper where we observed 
experimentally. As the tilt increases, we predict the on- 
set of stronger and stronger circulation. In this respect 
the predicted behavior the simulation and experiment are 
similar for small tilt angles. 

As the tilt increases further, we see in the simulations 
that the circulation cells move relative to the point on 
the hopper wall where we have observed the flow experi- 
mentally. This movement of the circulation cells puts the 
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FIG. 8: Ratio of azimuthal to radial velocity for a bare hopper 
as a function of hopper tilt angle e as imaged through the 
right (•) and left (o) windows of the test hopper, rotated 
to the frame where the untilted flow is entirely radial and 
mapped so that all tilt angles are positive. Lines indicate 
numerical predictions for von Mises (dotted) and Matsuoka- 
Nakai (dashed) plasticity models. 



observation area at a point where the rotational flow is 
much smaller and approaches zero. This decrease in the 
flow is not observed experimentally. 

We find that the values predicted using the Matsuoka- 
Nakai constitutive relation match our observations rea- 
sonably well for small angles, |e| < 0.1, while the gran- 
ular von Mises relation predicts much too little circula- 
tion. This suggests that the Matsuoka-Nakai criterion 
may more accurately describe the yielding of dense gran- 
ular flows in the regime we are studying. 

When the hopper wall friction is lowered (e.g. FigureO 
teflon lining) the ratio v$jv r is numerically predicted to 
be very small for both the Matsuoka-Nakai and von Mises 
yield criteria, as shown in Fig. [9l Both sets of predictions 
are beneath the threshold of detection in this regime and 
we cannot distinguish between them experimentally. We 
note, however, that we observe no trend for |e| < 0.2 in 
Fig. [7] 



VI. VELOCITY FLUCTUATIONS 

Although the numerical approach here describes the 
mean flow, we know that stress in hoppers is subject to 
large fluctuations about the mean Since Levy's flow 
rule assumes that the velocity in our hopper is directly 
related to the stress, it is interesting to examine any fluc- 
tuations in our velocity distributions [ItI ]. 

Previous work by Zhu and Yu using numerical Discrete 
Element Methods has studied the probability distribu- 
tion of velocities in a flat-bottomed cylindrical hopper as 
wall friction was varied from \i = 0.1 to 0.5 [13]. Al- 
though their cylindrical system had notable differences 
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FIG. 9: When wall friction is extremely low, /j, m — 0.05, the 
numerically-predicted ratio of azimuthal to radial velocity as 
a function of hopper tilt angle for both the Matsuoka-Nakai 
(o) and von Mises (□) criteria is four orders of magnitude 
smaller than the experimentally observed ratio at large tilt 
angles. 



from ours, including the presence of stagnant regions and 
plug flow, separate distributions for the radial component 
of flow were calculated, allowing rough comparison. In 
the DEM simulations it was found that the fluctuations 
decrease as wall friction is reduced. 

In Fig. [TU] we show the probability density distribu- 
tion of radial velocity in an untilted hopper for the four 
wall frictions we examined. We generate our velocity 
histograms by scaling several sets of data by the mean 
radial velocity for that set. Since the mean velocity var- 
ied for some experiments, for each value of wall friction 
we use the largest data set with consistent radial veloc- 
ities to generate our histograms. Within one data set, 
the velocity fluctuates sharply with time due to both our 
tracking technique which can produce frames in which no 
velocities are assigned and the large variability inherent 
to dense granular flows [2(J HH • The power spectrum for 
the time series does not indicate periodicity. 

As friction is increased from moderate values we ob- 
serve a decrease in the overall velocity fluctuations in 
agreement with the numerical predictions. The coeffi- 
cient of friction for Teflon is below the range studied 
previously, and in that case we observe that the fluctua- 
tions are suppressed relative to moderate wall frictions, 
indicating a non-monotonic dependence of velocity fluc- 
tuations on wall friction. All our observed distributions 
have skew, positive for moderate friction and negative for 
extreme frictions, but because we have studied velocities 
interpolated from many frames we do not see the wide 
distributions of Moka and Nott [22] . 

Finally, we examined the distribution of velocities as 
a function of tilt angle. For all wall linings, independent 
of the presence of secondary circulation, we observed an 
increase in fluctuations around the mean values as the 
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FIG. 10: The width of the distributions of probability density for radial velocity varies non-monotonically with wall friction, 
with the extremes having much narrower distributions than more moderate friction values. The distributions are roughly fit 
by Gaussians (clashed lines), with moderate values of friction having positive skewness and extreme values having negative 
skewness. 



tilt angle was increased. Figure [TT] shows the probability 
distributions for the radial and azimuthal components of 
the velocity as well as the ratio of those two components 
for untilted and highly tilted cases. For both tilted and 
untiltcd hoppers the means are similar, but the distri- 
bution grows much wider with tilt angle. Interestingly, 
although the distributions are wider for the highly tilted 
hopper, the distribution of v^/vr is of roughly the same 
width as that of v^, indicating that the fluctuations of 
the two components are correlated. This is associated 
with the fact that individual grains follow roughly fixed 
trajectories even if they are different from the mean flow. 
These persistent trajectories may indicate long-term cor- 
relation of particle contacts that might be expected for 
clusters of grains moving together, such as the "granular 
eddies" of Erta§ and Halsey [23[ . 



VII. CONCLUSION 

We have observed that real three-dimensional flows of 
granular matter are more complex than the idealized case 
of perfectly radial Jenike-like flow. In response to small 
perturbations (in the form of tilts), non-radial velocities 
arise leading to potentially large shifts in behavior. For 
small tilts we have found that the magnitude of circu- 
lation numerically predicted using the Matsuoka-Nakai 
constitutive relation is closer to observations than nu- 
meric predictions made with the more traditional, gen- 
eralized von Mises condition. At larger tilts we observe 
that the non-radial flow becomes even more pronounced, 
exceeding numerical predictions for all constitutive rela- 
tions. 

We have examined the influence of wall friction and dis- 
covered that extreme values, low or high, act to suppress 
secondary circulation. In the case of large perturbations 
and low friction, we see an abrupt and unpredicted onset 
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of secondary circulation in the opposite direction of the 
moderate friction case. 



Finally, we have observed that the fluctuations of in- 
dividual grain velocities grow as tilt is increased. We 
find this to be true independent of lining materials. We 
also note that a single grain may follow a trajectory that 
is quite different from the mean flow for long periods of 
time, which produces large, correlated fluctuations. Such 
particle-scale fluctuations are not captured by continuum 
models, but do not seem to affect the mean behavior. 
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FIG. 11: Probability density distributions of a) the azimuthal 
(left) and radial (right) components of the velocity and b) 
^ for an untilted hopper (solid line) and a hopper at tilt 
e = 0.14 (dashed line). These results are for an Emery 
paper-lined hopper, although they are representative of other 
wall frictions. 
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